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Abstract 

The ground state energy per particle of a dilute, homogeneous, 
two-dimensional Bose gas, in the thermodynamic limit is shown rig- 
orously to be Eq/N = (2irh 2 p/m)\ \n(pa 2 )\~ l , to leading order, with a 
relative error at most O (| ln^a 2 )!" 1 / 5 ) . Here N is the number of par- 
ticles, p = N/V is the particle density and a is the scattering length 
of the two-body potential. We assume that the two-body potential 
is short range and nonnegative. The amusing feature of this result is 
that, in contrast to the three-dimensional case, the energy, Eq is not 
simply N(N — l)/2 times the energy of two particles in a large box of 
volume (area, really) V. It is much larger. 



1 Introduction 

An ancient problem, going back to the 1950 's, is the calculation of the ground 
state energy of a dilute Bose gas in the thermodynamic limit. The particles 

1 Work partially supported by U.S. National Science Foundation grant PHY 98-20650. 
© 1999 by the authors. This paper may be reproduced, in its entirety, for non-commercial 
purposes. 
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are assumed to interact only with a two-body potential and are enclosed in 
a box of side length L. A formula was derived for the energy Eq(N,L) in 
three dimensions for a two-body potential v with scattering length a (see 
Appendix) and fixed particle density p = N/V, (N = particle number and 
V = volume =L 3 in three dimensions). In the thermodynamic limit, the 
energy/particle is 



e(p)= Urn E {N,p" 1/3 N 1/3 )/N -Anfipa (1.1) 

N^oo 

to lowest order in p. Here, p = h 2 /2m with m the mass of a particle. 

Our goal here is to derive the analogous low density formula for a two- 
dimensional Bose gas. 

There were several approaches in the 50's and 60's to the derivation of the 
three-dimensional formula fll.lj) , but none of them were rigorous. Recently we 
were able to give a rigorous derivation of fll.lj ) and we refer the reader to [|TJ for 
a physically motivated discussion of the essential difficulty in proving ( | . 1|) , 
which, basically, is the fact that at low density the mean interparticle spacing 
is much smaller than the mean de Broglie wavelength of the particles. Thus, 
Bose particles cannot be thought of as localized. Furthermore, in |J, we 
explain rather carefully why the usual expression 'perturbation theory' is not 
appropriate for (JL.1|) — especially in the hard core case. Indeed, Bogolubov's 
1947 'perturbation theory' || yields an estimate, which is incorrect for the 
low density limit: 



<P) -\p\ v. (1.2) 

It was only with a leap of faith that Bogliubov and Landau recognized that 
J v is the first Born approximation to 87ipa and thus were able to derive ( |1.1Q . 
Obviously this cannot be called perturbation theory. Moreover, depending 
on the nature of v, it is sometimes the potential energy and sometimes the 
kinetic energy that is the dominating quantity; for example, in the hard core 
case the kinetic energy is the perturbation, rather than the potential energy, 
as the Bogolubov method assumes. 

The two-dimensional theory, in contrast, began to receive attention only 
much later. The first derivation of the correct asymptotic formula was, to 
our knowledge, done by Schick || for a gas of hard discs: 
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e(p) ~47r/ip|ln(pa 2 )|" 1 . (1.3) 

This was accomplished by an infinite summation of 'perturbation series' 
diagrams. Subsequently, a corrected modification of |3|] was given in Q. 
Positive temperature extensions were given in || and in ||. All this work 
involved an analysis in momentum space — as was the case for ( |1 . 1| ) , with the 
exception of a method due to one of us that works directly in configuration 
space J/J. Ovchinnikov |J derived ( |1.3| ) by using, basically, the method in 
0. Again, these derivations require several unproven assumptions and are 
not rigorous. 

One of the intriguing facts about (|1.3|) is that the energy for N particles is 
not equal to N(N — l)/2 times the energy for two particles in the low density 
limit — as is the case in three dimensions. The latter quantity, E (2,L), 
is, asymptotically for large L, equal to 87r/iL~ 2 [ln(L 2 /a 2 )] \ Thus, if the 
N(N — l)/2 rule were to apply, ( |1.3| ) would have to be replaced by the much 
smaller quantity 47rpp [ln(L 2 /a 2 )] . In other words, L, which tends to oo in 
the thermodynamic limit, has to be replaced by the mean particle separation, 
p- 1 / 2 j n the logarithmic factor. Various poetic formulations of this curious 
fact have been given, but the fact remains that the non-linearity is something 
that does not occur in more than two- dimensions and its precise nature is 
hardly obvious, physically. This anomaly is the main reason that the present 
investigation is not a trivial extension of Jl| . 

We will prove ( |1.3| ) for nonnegative, finite range two-body potentials by 
finding upper and lower bounds of the correct form. The restriction to finite 
range can be relaxed somewhat, as was done in ||, but the restriction to 
nonnegative v cannot be removed in the current state of our methodology. 
The upper bounds will have relative remainder terms 0(| ln(pa 2 ) while the 
lower bound will have remainder 0(| \n.{pa 2 )\~ 1 ^). It is claimed in M that the 
relative error for a hard core gas is negative and 0(ln | ln(pa 2 )|| ln(pa 2 )| _1 ), 
which is consistent with our bounds. 

In the next section we shall give the upper bound (following Dyson's anal- 
ysis for the three-dimensional hard core gas |10|)- Then we shall recall our 
method in for the lower bound and show how it has to be modified. An 
important point concerns the definition of the scattering length in two di- 
mensions (which will be discussed in detail in Appendix A) and how 'Dyson's 
Lemma' |T0| , |l| has to be modified accordingly (Appendix B). 
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An obvious extension of the present work is the case of 2D bosons in a trap 
and this will be the subject of a forthcoming paper. Just as the passage from 
3D to 2D for the homogeneous case presents some non-trivial issues that have 
to be resolved, so the correct generalization of the Gross-Pitaevskii equation 
| ill to the 2D dilute trapped gas presents some additional complications. 

We thank P. Kevrekidis for drawing our attention to this problem. 



2 Upper Bound for the Ground State Energy 

We begin with the well known definition of the Hamiltonian under discussion: 



N 



ffM = -Ai£v? + $>(!*« -sjl), (2.1) 

i=l i<j 

We assume that v(r) > and v(r) = if r > Rq, for some Rq < oo. The 
Hamiltonian ( |2.1| ) acts on totally symmetric, square integrable wave functions 
of (xi, . . . ,xn) with Xi G M 2 . Its ground state energy in a box (rectangle, 
actually) of side length L is 

E (N, L) = inf (2.2) 

where the infimum is over all wave functions ^ satisfying appropriate con- 
ditions on the boundary of the box. For the upper bound it is natural to 
use Dirichlet boundary conditions, which gives the largest energy, but for the 
actual calculations it is more convenient to use periodic boundary conditions 
and a periodic extension of the interaction potential. This can only raise the 
energy since v > 0. Localization of the wave functions on the length scale L 
to obtain Dirichlet boundary conditions costs an energy ~ ((const.) L~ 2 per 
particle, so in the thermodynamic limit our upper bound is also a valid upper 
bound for Dirichlet boundary conditions. For the lower bound, on the other 
hand, we shall use Neumann boundary conditions, which yields the smallest 
energy. 



Following |nj we make a variational ansatz for \1/ of the following form: 

N 

*(x 1 ,...,x N ) = Y[f(t i (x 1 ,...,x i )) (2.3) 

i=2 
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where £, = min{|xi — Xj\, 1 < j < i — 1} is the distance of x« to its nearest 
neighbor among the points 2i, . . . , Xj_i and / is a nondecreasing function of 
t > with values between zero and 1. 

We wish to calculate (\I>, if (Ar) \[>)/(\[>, \[>). Dyson carried out this 
calculation for the hard core case, namely when f(r) = for r < the core 
radius. His formula has been generalized in || in two directions: One is the 
inclusion of an external potential (which we do not need here) and the other 
(which we do need) is the extension to a non-hard core potential v. We refer 
to §J for details. The result involves the following three integrals 

POO 

I = 2ir / (1 - f(r) 2 )rdr (2.4) 
Jo 

POO 

J = 2vr / (\f(r)\ 2 + lv(r)\f(r)\ 2 )rdr (2.5) 
Jo 

POO 

K = 2ir / f(r)f'(r)rdr (2.6) 
Jo 

In terms of these integrals the bound on the energy is 

(vP, VP) < N (jd- + Ij^y) ■ (2-7) 

The form of this bound is the same as in fllQf . Compared to Eq. (3.29) in 
there is a factor (1 — pl)^ 1 in the first term in place of (1 — pl)~ 2 - This can 
be traced to the use of the Cauchy-Schwarz inequality in Eq. (3.19) in 
which is not necessary in the case of the homogeneous system treated here. 

The next step is to make a choice for /, and this will involve the scattering 
length a and a variational parameter b. First, we have to define the scattering 
length. 

Consider the Schrodinger equation 

-fiAfa + \v<p = 0. (2.8) 

We do not require <po to be bounded. As shown in Appendix A, up to an 
overall factor there is a unique, nonnegative, spherically symmetric <po(x) = 
fo(\x\) that satisfies ( |2.8| ) provided the Schrodinger operator — pA + \v{r) in 
L 2 (R 2 ) has no bound states. For r > Rq, f necessarily has the form (since 
0o is a harmonic function outside the range of v) 



fo{ r ) — (const.) ln(r/a). 



(2.9) 
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The length a is called the scattering length. Note that it depends on both fi 
and on v. In the case that v is nonnegative, /o is necessarily a monotonically 
increasing function of r. 

We now define our variational / to be 

f(T) = f/.M/A(6) for0<r<6, 
1 1 for r > b. 

with some b > R Q > a to be chosen in an optimal way. By Appendix A we 
have that / satisfies /' > and < / < 1, for all b. Moreover, 

f( r ) > | ln ( r / a )/ ln ( 6 / a ) ioi a<r<b, 
~ 1 for r < a. 

Using this information one computes 

1 * - 2 + 2 ^X 1 -^)*^ (1+0(W6/a)1 " 1)) 

(2.12) 

J = ^^rW!-^ (2.13) 



K = ir j (f(r) 2 )'r dr = nb — tx J f(r) 2 dr 

£ rt -»jf^*=h&( 1 + °(WV.)]-'). (2-14) 

Inserted in (|2.7| ) this leads to the upper bound 

The minimum over b of the leading term is obtained for b = {2ixp)~ l l 2 . 
Inserting this in (|2.15| ) we thus obtain 

Theorem 2.1 (Upper bound). The ground state energy with periodic 
boundary conditions satisfies 

E (N, L)/N < j^y { (1 + 0(| Mpa 2 )\-i) . (2.16) 

Dirichlet boundary conditions may introduce an additional relative error, 
but as already noted it is at most AEq/N oc L~ 2 . 
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3 Lower Bound to the Ground State Energy 

The method of ||] for obtaining a lower bound to Eq(N,L) involves the 
following steps: 

1 A generalization of a lemma due to Dyson that allows the replace- 
ment of the interaction potential v by a 'soft' potential U at the cost 
of sacrificing kinetic energy. 

2 Division of the large box of side length L into small boxes of side 
length £, which is kept fixed as L — > oo, and a corresponding lowering 
of the energy by the use of Neumann boundary conditions on each 
box. It is necessary to minimize the total energy over all distributions 
of the particles among the small boxes; this is accomplished with the 
aid of the superadditivity of the ground state energy in each box (i.e., 
£ (iVi + N 2 ,L)> E (N U L) + E (N 2 , L), which follows from v > 0). 

3 The use of a rigorous version of first order perturbation theory, known 
as Temple's inequality [12|], to estimate from below the energy with the 
new potential U in the small boxes. 

We follow the same strategy here, but there are several modifications to 
be made. The two dimensional version of the generalized Dyson Lemma is 
as follows. 

Lemma 3.1. Let v(r) > and v(r) = for r > Rq. Let U(r) > be any 
function satisfying 

U(r)\n(r/a)rdr < 1 and U(r) = for r < Rq. (3.1) 

Let B C M 2 be star-shaped with respect to 0. Then, for all functions <p £ 
H\B), 



/i|V0(x)r + -v(r)\<p(x)\ z d 2 x>Li / U(r)\(j)(x)\ 2 d l x. (3.2) 

A domain B is star-shaped with respect to a point p if the line segment 
\p,x] C B whenever x G B. A convex domain is star-shaped with respect to 
any point in it (and conversely). The three-dimensional version of the lemma 
replaces ( [3.1|) with J °° U{r)r 2 dr < a. 
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The proof is given in Appendix B. 

As in [H], Lemma [3J] can be used to bound the many body Hamiltonian 
Hi 1 *) from below, as follows: 



Corollary 3.1. For any U as in Lemma \3. i| and any < e < 1 

H (N) > £T (iv) + (1 _ (3.3) 

W(xi, . . . , xn) = U I min \x{ — xA. I . (3.4) 

For [/ we choose the following functions, parameterized by R > Rq: 

U K (r) = h R ^ f ° r R ° <r<R (3.5) 
I otherwise 

with v(R) chosen so that 

rR 

/ Ur\t) ln(r/a)r dr = 1 (3.6) 



for all i? > i?o, i.e., 

= / In(r/o)r dr = \ {R 2 (\n{R 2 /a 2 ) - l) - R 2 (ln(R 2 /a 2 ) - l) } . 

JRo 

(3.7) 



The nearest neighbor interaction ( |3.4|) corresponding to Ur will be denoted 
W R . 

As in Jl| we shall need estimates on the expectation value, (Wr)o, of Wr 
in the ground state of eT( n ' of ( |3.3j ) with Neumann boundary conditions. 
This is just the average value of Wr in a hypercube in M. 2N . Besides the 
normalization factor v(R), the computation involves the volume (area) of 
the support of Ur, which is 

A{R) = tt{R 2 -R 2 q ). (3.8) 

In contrast to the three-dimensional situation the normalization factor 
u(R) is not just a constant (R independent) multiple of A(R); the factor 
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ln(r/a) in (|3.1|) accounts for the more complicated expressions in the two- 
dimensional case. Taking into account that Ur is proportional to the charac- 
teristic function of a disc of radius R with a hole of radius i?o, the following 
inequalities for n particles in a box of side length I are obtained by the same 
geometric reasoning as in [p]]: 

(Wr) > ^ (1 - f ) 2 [1 - (1 - Q) (n - 1} ] (3.9) 

(Wr)o < j^r [1 - (1 - Q) in - 1] } (3.10) 

with 

Q = A{R)/f (3.11) 

being the relative volume occupied by the support of the potential Ur. Since 
U R = v(R)~ 1 Ur we also have 

(W 2 R ) < -j^AWrW (3.12) 
As in [|] we estimate [1 - (1 - Q) (n ~ 1} ] by 

<„- i W >[i-<i > T *Lzi^ (3,3) 

This gives 

<*■>■ >- ^'Tr£m> w 

(w R )o < ^n^-Q- (3-15) 

From Temple's inequality (see H, |§) we obtain the estimate 

/ . - ^ i Ti- - I , K( w r)o ~~ ( W r)o) \ l0 ...... 

%^ > 1-ePso 1-— 77-1 . . : 3.16 



where 



S!°> = f (3.17) 
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is the energy of the lowest excited state of eT^ n \ This estimate is valid for 
Ei /n > (Wr)o, i.e., it is important that £ is not too big. 
Putting ( |3.14j ) and ( |3.16| ) together we obtain the estimate 



. „, n(n - 1) A(R) , . , 
E (n,£) > y & ' ^K(n) (3.18) 



with 



*(") = (!-')■ JTJ^Q 'y~ (ev(R)IP)-n(n-l)o ) ^ 

Note that Q depends on £ and R, and K depends on £, R and s besides n. 
We have here dropped the term (Wji)l in the numerator in ( ^.16| ), which is 
appropriate for the purpose of a lower bound. 

We note that K is monotonically decreasing in n, so for a given n we 
may replace K(n) by K(p) provided p > n. As explained in ^ convexity 
of n i— > n(n — 1) together with superadditivity of Eo(n,£) in n leads, for 
p = 4:p£ 2 , to an estimate for the energy of N particles in the large box when 
the side length L is an integer multiple of I: 

Eo(N, L)/N > (l - K(4pf ) (3.20) 

with p = N/L 3 . 

Let us now look at the conditions on the parameters e, and £ that have 
to be met in order to obtain a lower bound with the same leading term as 
the upper bound (|2.15| ). 

From ( |3.7|) we have 

W> = wm^) {1 -*™^ HR,B ° )) (3 - 21) 

We thus see that as long as a < R < p^ 1 ^ 2 the logarithmic factor in the 
denominator in (|3.20| ) has the right form for a lower bound. Moreover, for 
Temple's inequality the denominator in the second factor in ( p.!9| ) must be 
positive. With n = Ap£ 2 and u{R) > (const. )R 2 \n(R 2 /a 2 ) for R > R , this 
condition amounts to 



{const.)e\n(R 2 /a 2 )/£ 2 > p 2 £ 4 . 



(3.22) 
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The relative error terms in ( |3.20| ) that have to be 1 are 

, 1 R nR 2 Pt (W\ 

We now choose 

e~\\n(pa 2 )\-^, £ ~ p~ l ' 2 \ Hpa 2 )^ 10 , R ~ p^lMpa 2 )^ 10 

(3.24) 

Condition ( p.22|) is satisfied since the left side is > (const.) | ln(pa 2 )| 3//5 
and the right side is ~ | ln(pa 2 )| 2 / 5 . The first three error terms in (|3.23j ) are 
all of the same order, | ln(pa 2 )| -1 / 5 , the last is ~ | ln(pa 2 )|~ 1//5 (ln | ln(pa 2 )|) _1 . 
With these choices, (|3.20| ) thus leads to the following: 

Theorem 3.1 (Lower bound). For all N and L large enough such that 
L > (const. )p~ 1/2 | \n(pa 2 )\ l / w and N > (const.)] ln(pa 2 )| 1/5 with p = N/L 2 , 
the ground state energy with Neumann boundary condition satisfies 

E (N, L)/N > (1 - 0(| Hpa 2 )\-^)) . (3.25) 

In combination with the upper bound of Theorem 2.1 this also proves 
Theorem 3.2 (Energy at low density in the thermodynamic limit). 



lim — - = 1 (3.26) 

pa^o Anpp\ h^pa 2 )!" 1 

where eo(p) = liniAr^oo E (N, p~ l l 2 N l l 2 ) / N . This holds irrespective of bound- 
ary conditions. 

Remarks: 1. It follows from the the remark at the end of Appendix A 



that Theorem |3]2| is also valid for an infinite range potential v provided that 
v > and that for some R we have v(r)r dr < oo. 

2. As in || we could derive explicit bounds for the error term in 
flOp , but there is little reason to belabor this point. 
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A Appendix: Definition and Properties of 
Scattering Length 

In this appendix we shall define and derive the scattering length and some 
of its properties. The reader is referred to |T3| , especially chapters 9 and 11, 
for many of the concepts and facts we shall use here. While we are interested 
in two dimensions, much of the following is valid in all dimensions. 

We start with a potential \v(x) that depends only on the radius, r = \x\, 
with x G W 1 . For simplicity, we assume that v has finite range; this condition 
can easily be relaxed, but we shall not do so here, except for a remark at 
the end that shows how to extend the concepts to infinite range, nonnegative 
potentials. Thus, we assume that 

v(r) = for r > R . (A.l) 

We decompose v into its positive and negative parts, v = v + — t>_, with 
v + , i>_ > 0, and assume the following for t>_ only (with e > 0): 

'L 1 ^ 1 ) for n = l 
V- G { L 1+e (R 2 ) torn = 2 (A.2) 
L n / 2 (M n ) forn>3. 

In fact, v can even be a finite, spherically symmetric measure, e.g., a sum of 
delta functions. 

We also make the important assumption that \v(x) has no negative 
energy bound states in L 2 (M n ), which is to say we assume that for all 
cf) G iJ 1 (lR n ) (the space of L 2 functions with L 2 derivatives) 

/i|V0(x)| 2 + ^v(x)\4>(x)\ 2 d n x > (A.3) 



Theorem A.l. Let R > R and let B R C W 1 denote the ball {x : < \x\ < 
R} and Sr the sphere {x : \x\ = R}. For f G if 1 (i?^) we set 

S R [<P]= [ //|V0(x)| 2 + ^(x)|^)| 2 (A.4) 

J Br 1 
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Then, in the subclass of functions such that <p(x) = 1 for all x G Sr, there 
is a unique function (fio that minimizes 8r[<P]. This function is nonnegative 
and spherically symmetric, i.e, 

o (x) = / (N) (A.5) 

with a nonnegative function /o on the interval (0, R], and it satisfies the 
equation 

-}tA(f) (x) + ^v{x)(f) {x) = (A.6) 

in the sense of distributions on Br, with boundary condition fo(R) = 1. 
For R <r < R 

(r — a)/(R — a) for n = 1 

Mr) JT ym 'ir) = { ln(r/a)/HR/a) forn = 2 (A.7) 

[1 - ar 2 ~ n ) /(l - aR 2 - n ) for n > 3 

/or some number a called the scattering length. 
The minimum value of £R[<f)] is 

(2fx/(R-a) forn = l 

2nfi/\n(R/a) for n = 2 (A.8) 

27r rt / 2 /ia/[r(n/2)(l - aR 2 - n )} for n > 3. 



Remarks: 1. Given that the minimizer is spherically symmetric for 
every R, it is then easy to see that the R dependence is trivial. There is 
really one function, F , defined on all of the positive half axis, such that 
fo{ r ) — -^o( r )/-^o(-R)- That is why we did not bother to indicate the explicit 
dependence of fo on R. The reason is a simple one: If R > R, take the 
minimizer f for R and replace its values for r < R by fo(r)f (R), where fo 
is the minimizer for the Br problem. This substitution cannot increase 8^. 
Thus, by uniqueness, we must have that fo(r) = fo(r)fo(R) for r < R. 



From ( |A.7|) we then see that /o Symp (r) > for all r > Rq, which 



implies that a < Rq for n < 3 and a < Rq 2 for n > 3. 
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3. According to our definition (|A.7|) , a has the dimension of a length 
only when n < 3. 



4. The variational principle (|A.4|) , ( |A.8|) allows us to discuss the con- 



nection between the scattering length and J v. We recall Bogolubov's per- 
turbation theory 0, which says that to leading order in the density p, the 
energy per particle of a Bose gas is e (p) ~ 2irp J v, whereas the correct 
formula in two- dimensions is 4irpp\ ln(pa 2 )| _1 . The Bogolubov formula is 
an upper bound (for all p) since it is the expectation value of H^ N ' in the 
non-interacting ground state \P = 1. Thus, we must have \ J v > | ln(pa 2 )| _1 
when pa 2 <C 1, which suggests that 

" - sfe- (A ' 9) 



Indeed, the truth of ( |A.9| ) can be verified by using the function = 1 



as a trial function in ( |A.4| ). Then, using (|A.8|) , | J v > E = 2np/ \n(R/a) 
for all R > Rq, which proves (|A.9|) . As a — > 0, ( |A.9| ) becomes an equality, 
however, in the sense that ( L 2 v) \n(R /a) — > inp. 

In the same way, we can derive the inequality of Spruch and Rosenberg 



\TM for dimension 3 or more: 



• £ wk- (A ' 10) 



(Here, we take the limit R — »• oo in ( |A.8|) ). 
In one-dimension we obtain (with R = Rq) 



v > (A.ll) 
Ro — a 



Proof of Theorem \A.1\ : Given any G H 1 we can replace it by the square 
root of the spherical average of \<fi\ 2 - This preserves the boundary condition 
at \x\ = R, while the v term in ( |A.4p is unchanged. It also lowers the gradient 
term in (|A.4|) because the map p i— > J (Vy 7 /)) 2 is convex [13|. Indeed, there is 



a strict decrease unless <p is already spherically symmetric and nonnegative. 

Thus, without loss of generality, we may consider only nonnegative, spher- 
ically symmetric functions. We may also assume that in the annular region 
A = {x : Rq < \x\ < R} there is some a such that ( |A.7| ) is true because 



these are the only spherically symmetric, harmonic functions in A. If we 



LY 4/02/00 



15 



substitute for the harmonic function in A that agrees with at \x\ = Rq 
and \x\ = 1 we will lower Sr unless is already harmonic in A. (We allow 
the possibility a = for n < 2, meaning that = constant.) 

Next, we note that £r[0] is bounded below. If it were not bounded then 
(with R fixed) we could find a sequence J such that Sr^) — > — oo. However, 
if ft, is a smooth function on M + with h(r) = 1 for r < R+l and /i(r) = for 
r > 2i? + l then the function <j)i(x) = 0- ? (x) for \x\ < R and <p^{x) = h(\x\) for 
|x| > R is a legitimate variational function for the L 2 (M n ) problem in ( |A.3[ ). 
It is easy to see that £r[0 j ] < £r[&} + (const) R n ~ 2 , and this contradicts 
flAll ) (recall that is fixed). 

Now we take a minimizing sequence J for £r and corresponding 0-? as 
above. By the assumptions on u_ we can see that the kinetic energy T- 7 = 
J" |V0 J | 2 and f |0 J | 2 are bounded. We can then find a subsequence of the J ' 
that converges weakly in H 1 to some spherically symmetric 4>o(x) = fo(\x\). 
Correspondingly, J (x) converges weakly in H l {B R ) to <f> Q {x) = fo(\x\). The 
important point is that the term -J t>_|0 J | 2 is weakly continuous while the 
term f t> + |0 J | 2 is weakly lower continuous |L3||. We also note that fo(R) = 1 
since the functions $ are identically equal to 1 for R < \x\ < R + l and the 
limit 0o is continuous away from the origin since it is spherically symmetric 
and in H 1 . 

Thus, the limit function 0o is a minimizer for S[4>] under the condition 
= 1 on Sr. Since it is a minimizer, it must be harmonic in A, so (|A.7|) 
is true. Eq. ( |A.6| ) is standard and is obtained by replacing O by 0o + Sip, 
where ip is any infinitely differentiable function that is zero for \x\ > R. The 
first variation in 5 gives (|A.6| ). 

Eq. ( |A.8| ) is obtained by using integration by parts to compute £r[4>q\. 

The uniqueness of the minimizer can be proved in two ways. One way 
is to note that if 0o 7^ tpo are tw° minimizers then, by the convexity noted 
above, Sr[\/ 0q + ipo) < ^i?[0o] +£r{'4 } o)- The second way is to notice that all 
minimizers satisfy ( |A.6| ), which is a linear, ordinary differential equation for 
/o on (0, R) since all minimizers are spherically symmetric, as we noted. But 
the solution of such equations, given the value at the end points, is unique. 
□ 

We thus see that if the Schrodinger operator on MJ 1 with potential \v(x) 
has no negative energy bound state then the scattering length in ( |A.7| ) is well 
defined by a variational principle. Our next task is to find some properties 
of the minimizer 0q. For this purpose we shall henceforth assume that v is 
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nonnegative, which guarantees (|A.3|), of course. 



Lemma A.l. If v is nonnegative then for all < r < R the minimizer 
M x ) = fo(\x\) satisfies 
A) 

fair) > / aSymP W, (A.12) 

where f^ ymp is given in ( \A.% ) 

B) fo{ r ) is a monotonically nondecreasing function ofr. 

C) If v(r) > v(r) > for all r then the corresponding minimizers satisfy 
fo{ r ) < fo( r ) f or oil r < R. Hence, a > a > 0. 

Proof. Let us define f^ ymp (r) for all < r < oo by ( |A.7| ), and let us extend 
fo(r) to all < r < oo by setting /o( r ) — fo hymp ( r ) when r > R. 

To prove A) Note that — A0o = —\v(fio, which implies that 4>o is sub- 
harmonic (we use v > and 0o > 0, by Theorem [Al|) . Set h £ (r) = 
/o( r ) — (1 + s)fo t ' ymp (r) with e > and small. Obviously, x t— > h e (\x\) is 
subharmonic on the open set {x : < \x\ < oo} because /o Symp (M) i s har- 
monic there. Clearly, h £ — > — oo as r — > oo and h e (R) = —e. Suppose that 
( A.12 ) is false at some radius p < R and that ho(p) = — c < 0. In the annulus 



p < r < oo, h £ (r) has its maximum on the boundary, i.e., either at p or at 
oo (since /i(|x|) is subharmonic in x ). By choosing e sufficiently small and 
positive we can have that h e (p) < —2e and this contradicts the fact that the 
maximum (which is at least — e ) is on the boundary. 

B) is proved by noting (by subharmonicity again) that the maximum of 
fo in (0, r) occurs on the boundary, i.e., /o( r ) > fo( r ') for any r' < r. 

C) is proved by studying the function g = f — f . Since f and f Q are 
continuous, the falsity of C) implies the existence some open subset, fl C Br 
on which g(|o;|) > 0. On Q we have that g(|x|) is subharmonic (because 
v fo > vfo). Hence, its maximum occurs on the boundary, but g = there. 
This contradicts g(\x\) > on fl. 

□ 

Remark about infinite range potentials: If v(r) is infinite range and non- 
negative it is easy to extend the definition of the scattering length under the 
assumptions: 

1) v(r) > for all r and 

2) For some R\ we have fft 1 'v(r)r n ~ 1 dr < oo. 
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If we cut off the potential at some point Ro > R± (i.e., set v(r) = 
for r > Rq) then the scattering length is well defined but it will depend 
on R , of course. Denote it by a(R ). By part C of Lemma ( |A.1| ), a(R ) 
is an increasing function of R . However, the bounds ( |A.9| ) and (|A.10| ) and 
assumption 2) above guarantee that a(R ) is bounded above. (More precisely, 
we need a simple modification of (|A.9| ) and ( |A.10| ) to the potential v(r) = oo 
for r < Ri and v(r) = v(r) for r > R\. This is accomplished by replacing 
the 'trial function' f(x) = 1 by a smooth radial function that equals for 
r < R\ and equals 1 for r > R 2 for some R 2 > Ri-) Thus, a is well defined 
by 

a = lim a{Ro) . (A.13) 



B Appendix: Proof of Dyson's Lemma g7T] in 
Two Dimensions 

Proof. In polar coordinates, r,9, one has |V0| 2 > \d<p/dr\ 2 . Therefore, it 
suffices to prove that for each angle 6 G [0,2-71"), and with 4>(r,9) denoted 
simply by /(r), 

R(6) 1 f R(8) 

fi\df(r)/dr\ 2 + -v{r)\f(r)\ 2 rdr > fi U{r)\f(r)\ 2 rdr, (B.l) 

2 Jo 

where R{0) denotes the distance of the origin to the boundary of B along the 
ray 9. 

If R(6) < -Ro then ( |B.1| ) is trivial because the right side is zero while the 
left side is evidently nonnegative. (Here, v > is used.) 

If R{9) > Rq for some given value of 9, consider the disc V(9) = {x G 
M? : < \x\ < R{9)} centered at the origin in R 2 and of radius R(9). Our 
function / defines a spherically symmetric function, x h- > f(\x\) on T>(9), and 
( |B.1| ) is equivalent to 

/^|V/(M)| 2 + kr)|/(M)| 2 (B.2) 

V(6) A 



Now choose some R G (Ro, R(9)) and note that the left side of (|B.2|) is 
not smaller than the same quantity with T>(9) replaced by the smaller disc 
= {a; G M 2 : < \x\ < R}. (Again, v > is used.) According to 
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Appendix A, Theorem [A.l| , eq. ( |A.8| ), and linearity in \f\ 2 , this integral over 
V R is at least E(R)\f(R)\ 2 . Hence, for every R < R < R(8), 



The proof is completed by noting that E(R) = 27T/^/ ln(R/a), by multiplying 
both sides of ( |B.3| ) by U(R)R\ia(R/a) and, finally, integrating with respect 
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